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Superscaling of the quasielastic cross section in charged current neutrino-nucleus reactions at 
energies of a few GeV is investigated within the framework of the relativistic impulse approximation. 
Several approaches are used to describe final state interactions and comparisons are made with 
the plane wave approximation. Superscaling is very successful in all cases. The scaling function 
obtained using a relativistic mean field for the final states shows an asymmetric shape with a long tail 
extending towards positive values of the scaling variable, in excellent agreement with the behavior 
presented by the experimental scaling function. 



PACS numbers: 25.30.Pt; 13.15.+g; 24.10.Jv 



In the context of inclusive quasielastic (QE) electron 
scattering at intermediate to high energies, the concepts 
of scaling p| and superscaling Q have been explored in 
previous work 0, Q , where an exhaustive analysis of the 
(e, e') world data demonstrated the quality of the scal- 
ing behavior. Scaling of the first kind (no dependence on 
the momentum transfer) is reasonably well respected at 
excitation energies below the QE peak, whereas scaling 
of second kind (no dependence on the nuclear species) is 
excellent in the same region. The simultaneous occur- 
rence of both kinds of scaling is called superscaling. At 
energies above the QE peak both scaling of the first and, 
to a lesser extent, of the second kind are shown to be vi- 
olated because of important contributions introduced by 
effects beyond the impulse approximation, namely, in- 
elastic scattering p| together with correlations and me- 
son exchange currents in both the lp-lh and 2p-2h sec- 
tors |a,|7|. 

The scaling analysis of (e, e') data has recently been 
extended through the QE peak into the A region 
Of relevance to the present work we note that the high- 
energy inclusive electron scattering cross section is well 
represented up to the A peak using the scaling ideas, 
importantly, with an asymmetric QE scaling function. In 
that study the scaling approach was also used to predict 
nuclear (y, fi) cross sections, based on the assumption of 
a universal scaling function, valid for both electron and 
neutrino scattering at corresponding kinematics. 

In this letter we investigate the QE scaling proper- 
ties of charged-current (CC) neutrino-nucleus scattering 
within the context of the Relativistic Impulse Approxi- 
mation (RIA). After verifying that various RIA models 
do superscale, we compare the associated scaling func- 
tions with the (e, e') phenomenological one referred to 
above. This allows a check on the consistency of the uni- 



versality assumption and on the capabilities of different 
models to yield the required properties of the experimen- 
tal scaling function, specifically, its asymmetric form. 

Here we follow the general procedure of scaling and 
superscaling studies, namely we first construct inclusive 
cross sections within a model and then obtain scaling 
functions by dividing them by the relevant single-nucleon 
cross sections weighted by the corresponding proton and 
neutron numbers H,0,l3- The scaling function is plotted 
against the scaling variable ip{q,Lo), with q and ui the 
momentum and energy transferred in the process, and 
its scaling properties analyzed. 

Within the RIA framework CC neutrino-nucleus QE 
scattering is described by assuming that at the v-p, ver- 
tex one vector boson is exchanged with the nucleus, in- 
teracting with only one nucleon, which is then emitted 
while the remaining (A-l) nucleons in the target remain 
as spectators. The nuclear current operator is thus taken 
to be the sum of single-nucleon currents, for which we em- 
ploy the usual relativistic free nucleon expressions |8L flO| - 
The RIA approach has been extensively and successfully 
applied in investigations of exclusive electron scattering 
reactions [llf . Further details on the model have been 
presented in [13, [H G3 . 

We describe the bound nucleon states as self-consistent 
Dirac-Hartree solutions, derived within a Relativistic 
Mean Field (RMF) approach using a Lagrangian con- 
taining er, ui and p mesons 14j. For the description of 
the outgoing nucleon states we consider several different 
approaches. In one we use plane- wave spinors (thus no 
Final-State Interactions (FSI)), corresponding to the Rel- 
ativistic Plane- Wave Impulse Approximation (RPWIA). 
However, comparisons with data require a more realis- 
tic description of the final nucleon state, which should 
include the effects due to FSI. This is accomplished by 
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using solutions of a Dirac equation containing relativistic 
potentials. This constitutes the Relativistic Distorted- 
Wave Impulse Approximation (RDWIA). 

The use of complex relativistic optical potentials fit- 
ted to elastic proton scattering data has proven to be 
successful in describing exclusive (e, e'p) scattering reac- 
tions 11]. The absorption produced by the imaginary 
term represents the loss of flux into inelastic channels. 
For inclusive processes such as (e, e') and [y, fi), where a 
selection of the exclusive single-nucleon knockout chan- 
nel cannot be made, the contribution from these inelastic 
channels should be retained. Ignoring them would lead 
to an underestimation of the inclusive cross section |13| . 
A simple way of obtaining the right inclusive strength 
within the RIA is to use purely real potentials. We 
consider two choices for the real part. The first uses 
the phenomenological relativistic optical potential from 
the energy-dependent, A- independent parameterizations 
(EDAIC, EDAIO, EDAICa) derived by Clark et al. [H|, 
but with their imaginary parts set to zero. The sec- 
ond approach employs distorted waves obtained with the 
same relativistic mean field used to describe the initial 
bound nucleon states. We refer to these two FSI descrip- 
tions as real Relativistic Optical Potential (rROP) and 
Relativistic Mean Field (RMF), respectively. Dispersion 
relation and Green function techniques have been used 
to derive rigorously the potentials for inclusive scatter- 
ing (see refs. 0, [2, ^3) leading to results which are 
within a few percent of those obtained in the IA with ei- 
ther the rROP [l8| or the mean field [l6|. Note that the 
RMF model (1) is known to work quite well for inclusive 
QE (e, e') scattering (which is verified here; see below) 
and (2) is constructed in a way that fulfills the dispersion 
relation [l6| and maintains the continuity equation. 
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FIG. 1: Quasielastic differential cross section do / dE^dQ.^ 
versus the muon kinetic energy T M for the reaction (^ M ,/i _ ) 
on 12 C (left) and ie O (right). The incident neutrino energy 
is e v = 1 GeV and the muon scattering angle is (9 M = 45°. 
In each panel we present results for RPWIA (solid), rROP 
(dashed) and RMF (dot-dashed). The cross section for the 
RFG model is also presented for reference (dotted). 



Our results for CC neutrino-nucleus QE scattering are 
presented in Fig. ^ where we show the differential cross 
section (da/dE^dfl^) as a function of the outgoing muon 



kinetic energy for 12 C and 16 O. For reference we also in- 
clude the results obtained within the Relativistic Fermi 
Gas (RFG) (dotted) with Fermi momenta kp = 228 
MeV/c for 12 C 4j and k F = 216 MeV/c for 16 [ig|. The 
RFG curves contain a phenomenological energy shift 
yielding RFG cross sections which are very similar to the 
RPWIA results. 

The mean field dynamics in the initial and final nuclear 
states lead to cross sections having tails that extend both 
below and above the kinematical region where the RFG 
is defined. With FSI included we observe a reduction 
of the cross section, particularly in the case of the RMF 
potential where it is seen to be about 20% for both nu- 
clei in the region close to the maximum. Notice also a 
slight displacement in the maximum of the cross section 
in the cases of the two RIA-FSI models. However, the 
most striking feature is the long tail displayed by the 
RMF cross section for small muon kinetic energies. This 
corresponds to transferred energies above the QE peak, 
i.e., positive values of the scaling variable. These FSI ef- 
fects lead to a clear asymmetry in the RMF cross section, 
in contrast to the RPWIA and rROP results. The dis- 
crepancy between the two FSI approaches is linked to the 
different behaviors of the two potentials: for high nucleon 
kinetic energies (small muon energies) the scalar and vec- 
tor energy-dependent potentials of the rROP model are 
significantly reduced with respect to the RMF [ll|. The 
latter approach has strong scalar and vector potentials 
that are delicately balanced and that shift the strength 
towards higher energies. This explains why the rROP 
cross section is closer to the RPWIA case, and moreover, 
why the main difference between the two models is ob- 
served in the region of lower muon energies. Note that 
the asymmetric broadening of the RMF cross section is 
similar to that observed in non-relativistic models of the 
FSI, described by some as medium modifications of the 
p-h propagator jl9| . 

Let us now study the superscaling properties. We 
present results for the scaling function f(ip')i which is ob- 
tained by dividing the calculated differential cross section 
of Fig. n by the single-nucleon cross section as given in 
Eqs. (45,52,86-94) of 8]. This function is plotted against 
the shifted QE scaling variable ip' defined as 
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where A' = (a; — E s hift) /2itin , k = q/2m,N, t' = k 2 — A 
and £f = \/l + {kp/mji) 2 - 1. The energy shift E shift 
has been taken from [j]. Results correspond to fixed 
scattering muon angle 9^ = 45°, although similar scaling 
functions are obtained with other values. 

Scaling of the first kind is explored in Fig. [21 where we 
present f(i/j') for 12 C at three different values of the in- 
cident energy. Each panel in the figure corresponds to a 
different description of the FSI. As one can see, the scal- 
ing function for the RPWIA and rROP models shows 
a very mild dependence on the momentum transfer in 
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FIG. 2: Scaling function for three values of the incident neu- 
trino energy e„. Results correspond to 12 C and # M = 45°. 
Top, middle and bottom panels refer to RPWIA, rROP and 
RMF models (see text for details). 



FIG. 3: Scaling function for three nuclei: 12 C, 16 and 



Results correspond to £„ = 1 GeV and # p 
and bottom panels as in previous figure. 
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both positive and negative ip' regions. In the case of the 
RMF, a slight shift occurs in the so-called "scaling re- 
gion" ip' < 0, whereas for ip' positive the model breaks 
scaling at roughly the 30% level in the energy region ex- 
plored here. This is not in conflict with the experimental 
(e, e') data that indeed leave room for some breaking of 
first-kind scaling in this region, due partly to A produc- 
tion and partly to other contributions, such as MEC and 
their associated correlations in the 2p-2h sector 0. It 
is striking that the RMF model, in spite of being based 
on the impulse approximation, leads to the same kind of 
behavior which is apparently not reproduced by uncorre- 
cted models in the impulse approximation. Importantly, 
the RMF scaling function exhibits a significant asymme- 
try, being larger for positive ip', which persists for all 
neutrino energies (actually increasing with e v ). 

Scaling of the second kind is studied in Fig. [21 where 
the scaling function evaluated for three nuclei is pre- 
sented using the three models. The values of the Fermi 
momentum used range from Uf = 216 MeV/c for 16 
to k F = 241 MeV/c for 40 Ca 4]. The initial bound 
states have been obtained using the parameters of the 
set NLSH 20] . Results with other parameterizations 
are similar and do not change the general conclusions. 
As observed, the differences introduced by changing nu- 
cleus are small. We may conclude that within the present 



model scaling of second kind is very successful. Indeed, 
this is just what is seen experimentally, at least for ip' < 
where scaling of the second kind is excellent 0, ■ 

Finally, in Fig. 0] we compare the model superscaling 
functions with the averaged QE phenomenological func- 
tion obtained from the analysis of (e, e') data [3. Wi l2ll| . 

First we observe the symmetric character of the RP- 
WIA and rROP results, which clearly differ from the ex- 
perimental function. On the contrary, the RMF curve 
displays a pronounced tail that extends toward positive 
values of ip' , following closely the asymmetric behavior 
of the data and yielding excellent agreement with the 
phenomenological scaling function. In the light of the 
analysis in |8| this immediately implies that the RMF 
approach yields excellent agreement with the experimen- 
tal (e, e') inclusive cross section. The asymmetric shape 
of the RMF result constitutes a basic difference not only 
from the other two models explored in this work, but also 
from other modeling presented in the literature, such as 
those in , where the long tail in the superscaling func- 
tion is absent. Indeed, the CDFM model for correlation 
effects presented in .22] is manifestly symmetrical around 
the QE peak. It should be remarked that all of the curves 
in Fig. 0] essentially satisfy the Coulomb sum rule, i.e., 
they integrate to unity. 

The results in Fig. ^demonstrate that the RMF model 
is successful at incorporating important dynamical effects 
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FIG. 4: Scaling function evaluated within the RPWIA (solid), 
rROP (dashed) and RMF (dot-dashed) approaches compared 
with the averaged experimental function together with a phe- 
nomenological parameterization of the data (dotted). 



the data has usually been ascribed to the role played by 
ingredients beyond the mean field, such as short-range 
correlations and two-body currents 7]. We show here 
that it can be explained within the RIA framework re- 
sorting only to one-body excitations, provided strong rel- 
ativistic potentials are included in the model. It would 
be interesting to see how a more elaborated formalism 
that includes multinucleon excitations as in |17| compares 
with superscaling data. 

In summary, we have shown that superscaling is ful- 
filled to high accuracy within the present relativistic im- 
pulse approximation in the QE region, and that this holds 
for the three different descriptions of FSI considered here. 
The asymmetric shape and the long tail at positive ip'- 
values observed in the experimental scaling function is 
reproduced only by the RMF model. This result rein- 
forces our confidence in the adequacy of descriptions of 
FSI effects for inclusive (e, e') and [y, fi) reactions when 
based on the RMF approach. 
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